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Abstract 

Th«  purpoaa  of  this  Invsstlgstlon  is  to  devslop  on  intsrsctivs, 
ossr  orlsntsd  software  package  for  direct  digital  control  systen  design. 
The  batch  aode  prograB*  FOBTBAC,  developed  by  Professor  Brian  Porter  of 
the  Itelverslty  of  Salford,  England,  Is  the  source  of  the  coiiq>utatlonal 
algorlthaw  used.  This  report  details  the  Input, output,  and  program 
sequence  control  software  developed  to  produce  an  efficient,  user 
oriented  Interactive  package. 

The  package  Is  very  forgiving  of  user  errors  and  gives  the  user 
collate  control  over  what  data  to  Input,  what  data  will  be  output,  and 
lAlcih  parts  of  the  program  are  executed.  The  package  Is  quite  useful 
In  the  design  process  for  discrete-time  time-optimal  control  systems, 
where  many  possible  control  parameter,  variations  must  be  examined. 
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TSTESAC  - AN  INTERACTIVE  SOFTWARE  PACKAGE 


FOR  DIRECT  DIGITAL  CONTROL  DESIOI 


I.  Introduction 

In  thn  past  few  years  the  digital  computer  has  become  a valuable 
! asset  In  the  control  system  design  process.  With  Increased  Interest 

. In  modem  control  theories,  the  computer  Is  fast  becoming  a necessity. 

A portion  of  a computer  program  for  control  system  design  Is  going  to 
be  the  coding  dealing  with  Input,  output,  and  program  sequence  control. 

Hany  computer  programs  In  classical  and  a few  In  modem  control 
theory  now  exist  (Ref  1,2, 3, 4,  and  5),  but  most  can  be  quite  cuaibersome 
for  the  user.  These  programs  are  written  In  terms  of  contputatlonal  and 
computer  resource  efficiency,  often  at  the  expense  of  user  efficiency. 
Dher  efficiency  Is  defined  as  the  capability  of  a user  to  get  exactly 
what  he  needs  from  a program  execution  (no  more  and  no  leas)  with  the 
amount  of  work  and  expertise  In  the  program  operation. 

Them  is  no  unique  solution  to  a given  control  problem,  and  thus 
an  Itaratlve  procedure  with  one  or  more  design  parametem  nay  be 
required.  This  can  beat  be  achieved  in  an  Interactive  environment 
lAere  the  user  can  make  selective  changes  In  the  design  parameters 
I . baaed  upon  Intermediate  outputs.  The  problem  then  Is  to  develop  an 

! Interactive,  uaer  oriented,  program  for  the  I-O  (input-output)  and 

1 

program  aequance  control  functlona  which  can  be  put  around  any  aet  of 
control  ayatem  dealgn  theory  algorlthma. 


Th«  Problem 


Professor  Brlsn  Porter  of  the  University  of  Salford,  England,  has 
been  developing  a control  system  design  theory  based  on  discrete  state 
variable  feedback  to  achieve  a time  optimal  response  at  any  given 
sampling  rate.  The  algorithms  for  this  design  process  have  been 
furnished  by  Professor  Porter  and  the  Air  Force  Flight  Dynamics  Lab 
In  the  computer  program  FORTRAC.  This  batch  mode  program  which  was 
written  for  the  University  of  Salford  computer  and  adapted  to  run  on 
the  Wrlght-Patterson  AFB  computer  system  along  with  the  SCOPE  operating 
system  can  be  cumbersome  for  the  user. 

The  problem  focused  on  In  this  Investigation  Is  to  develop  the 
Interactive  1-0  and  program  sequence  control  software  for  the  design 
algorithms  of  FORTRAC  which  will  give  a user  oriented  design  package 
that  %rlll  operate  efficiently  on  the  CDC-6600/Cyber-74  computer  system 
at  Wrlght-Patterson  AFB,  Ohio.  The  capability  of  this  package,  INTERAC, 
Is  then  demonstrated  by  synthesizing  control  laws  for  a fly-by-wlre 
system  for  the  longitudinal  axis  of  the  C-141  and  F-4E  aircraft. 


Scope  and  Obiectlves 

The  software  was  designed  to  meet  the  following  characteristics: 

(1)  Input  timing  and  format  will  be  under  user  control.  (2)  Program 
termination  will  not  be  caused  by  any  format  errors  In  an  Input  sequence. 
In  case  of  such  an  error,  the  user  will  be  notified  and  be  asked  to 
reenter  that  Item.  (3)  Output  content  and  format  will  be  user  defined. 
(4)  Only  user  desired  elements  of  the  program  need  be  executed  and  In 
user  defined  order.  (5)  Program  directions  will  be  printed  for  the 
beginning  user  or  they  can  be  suppressed  for  the  experienced  user. 
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^ i (6)  Additional  algorithms  can  be  added  to  the  package  with  a minimum 

of  additional  programming.  (7)  The  1-0  subroutines  can  be  used  directly 
in  any  program  requiring  a user  oriented,  interactive  environment. 

(8)  The  package  will  operate  within  the  60,000  word  field  length 
sivallable  to  the  interactive  users  on  the  CDC  computers  at  Wrlght- 


Patterson  AIB. 

The  control  laws  for  the  C-141  and  F-4E  fly-by-wlre  systems  are 
obtained  for  the  longitudinal  axis  and  at  only  one  flight  condition. 


I 

A 
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Development 

Ch^ter  2 gives  an  overview,  of  the  FORTRAC  package  with  specific 
eaphasls  on  the  algorithms  Incorporated  into  INTEIL4C.  The  specialized 
INTBBAC  software  is  discussed  in  chapter  3.  The  mathematical  models 
for  the  C-141  and  the  F-4E  are  developed  in  Chapter  4.  The  chapter 
concludes  with  a discussion  of  the  control  laws  obtained  for  the  two 
aircraft.  Chapter  5 presents  the  conclusions  and  recoomiendationa  for 
Inprovlng  the  software  package  and  for  continued  control  system 
development  for  the  C-141  and  F-4E. 

The  actual  FORTBAN  Extended  source  code  for  the  specialized 
IHTBIAC  aoftware  is  given  in  Appendix  A,  and  a user's  giilde  for  the 
IMTEBAC  package  is  contained  in  Appendix  B.  Appendix  C is  a 
programmer 'a  guide  for  using  the  specialized -IMTERAC  subroutines  in 
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I ftimlshttd  under  contract  to  the  Air  Force  Flight  Dynamics  Lab.  This 

I package  Is  a collection  of  42  subroutines  tied  together  by  a master 


program  written  for  the  particular  design  problem  attempted.  The  master 
program  furnished  to  the  Flight  Dynamics  Lab  will  design  a discrete 
control  law  for  a sixth  order  or  smaller  system,  design  an  observer, 
and  run  a time  simulation  of  the  system  for  the  resulting  controller. 

The  theory  Is  good  for  any  order  system,  but  a larger  dimensioned 
aaatar  program  would  be  required  for  higher  order  systems. 

Table  1 shows  the  subroutines  from  FORTRAC  grouped  Into  four 
ganaral  functional  categories.  Most  of  the  siibroutlnes  listed  under 
control,  law  synthesis  and  utility  are  included  In  INTERAC,  and  their 
specific  function  will  be  discussed  In  this  chapter.  A description 
of  the  other  subroutines  can  be  obtained  from  Reference  6. 

FOKIRAC  Is  designed  to  take  the  continuous  time  description  of  a 
linear  system  (Fig  1) , and  to  synthesize  a control  law  for  any  of  the 
following  classes  of  problems:  (1)  discrete-time  time-optimal 
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regulator,  (2)  discrete-time  time-optimal  disturbance  rejector,  and 
(3)  discrete-time  time-optimal  tracker. 
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Figure  1.  Continuous  Tlae  Syati 
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* Th«a«  names  had  to  be  shortened  from  the 
orlglnel  names  In  Reference  6 In  adapting 

) the  package  to  operate  on  the  G)C  computer  : 

syatam  at  Wrlght-Patterson  AFB. 


The  type  of  problem  to  be  solved  and  the  dimensions  of  the  input 
•yetam  are  handled  by  RPDAIA.  The  I-O  of  the  matrix  data  are  handled 
by  MAIBIO  and  HATPRIO  respectively.  These  functions  will  be  handled 
in  IHTEkAC  by  the  specialized  INTERAC  software. 

Sampled  Data  Transformation 

The  first  stsp  in  any  class  of  problem  is  to  transform  the 
coatlnuous  system  of  Figure  1»  described  by  Bq  (1),  to  a sampled-data 
system  description»  Eq  (2). 


mm 


x(t)  = A^t)  Bu(t)  ♦ Rdt) 
^t)  = C^t) 


(1) 





o 


o 


x(k*1)  = f(T)x(k)  * G(T)u(k)  ♦ ^(T)dk) 
^k)  = C^k) 


(2) 


All  dlscrsttt  functions  of  k «r«  l^llsd  to  b«  functions  of  kT. 

TWO  major  subroutines  are  used  in  this  operation.  The  first  Is  EIGSTR, 

mliieh  also  calls  EI6VAV.  The  eigenvalues  and  eigenvectors  of  the 
continuous  system  are  determined  with  the  eld  of  Mtrlx  manipulation 
routines  from  a comoion  user  mathesMtics  library.  The  description  of 
the  matrix  manipulation  routines  required  for  FOSTBAC  and  INTESAC  are 
given  later  In  this  chapter. 

The  aeeond  major  subroutine  used  Is  SAMPLE,  which  generates  the 
discrete  matrices  F,(T),  G(T),  and  M(T)  from  A,  B,  and  R.  The 
transformation  can  be  accooqtUshed  using  Eqs  (3) , (4) , and  (5) . 


f(T)  = e-^ 

(3) 

,T 

S<T)  ' ^ g dt 

(A) 

DRO-)  = /e^tgdt 

(5) 

In  gansral,  an  exact  solution  to  the  above  equations  cannot  be  easily 
dstmxmlnsd.  To  faelUtata  tha  solution  the  diagonal  matrix^  which 
satisfies  Eq  (6)  is  produced. 


A = UAU*'  (6) 

U is  ths  matrix  of  aiganvactors  <ri>tained  from  EIGSTR.  At  this 


point  the  siibxoutlne  SAMPLE  Is  limited  to  systems  which  hsve  non* 


reposted  elgenvsiues.  Multiple  eigenvalues  would  produce  off  diagonal 
taxBs  in  ^ » thus  this  simplification  could  not  be  used.  The  equations 


The  confutation  of  becomes  a set  of  scalar  problems  0^'  , 

where  the  *s  are  the  eigenvalues  of  A which  appear  on  the  diagonal 


B,  and  R are  taken  outside  the 


la  the  disturbance  rejection  and  tracking  problems,  the  next 
atop  is  to  augment  the  eyetem  with  discrete- cine  Integrators  shown 


j 
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th«  trackttx  is  2.(1'^)  • Th*  xationala  and  spacific  raquiramants 

for  tha  usa  of  Incagral  faadback  is  glvan  In  Rafaranca  7 (papars  7 and  9) ^ 
Tha  cosnand  Input  matrix,  ^ Is  astabllshad  for  tha  input 

coaund,  v(k),  in  tha  augmantad  system  aquation.  Tha  augmented  dlscrate 
system  is  than  rapresantad  by  Eq  (11). 


«hara 


AF  = 


f 0 

,CT  I. 


AG  = 


G 

L0. 


ADR 


r DRi 

— ADE  = 

rg- 

L^J 

L-deJ 

xXk)  = 

Li(k)J 


Control  Law  Synthasis 

Bm  control  law  to  be  synthaslsad  is  given  by  Eq  (13). 


u(k)  = Kx'(k) 


xXk-l)  = AFxtk)  ♦ AG  u(k)  * ADR  dk)  * APE  v(k)  (11) 


(12) 


(13) 


first  stop  ossa  subroutine  TIAHFAB  to  transfom  the  natrlcas 


AF  and  AG  Into  Brunovsky  controllable  canonical  form.  Thla  la  dona 
oalng  tha  flrat  method  of  Aplavlch  (Ref  8:124-*126).  The  utility 
aubroutinaa  SETQUP  through  STQDOWN  (Table  I)  are  uaad  to  achieve  the 
txanafomatlon  and  to  generate  the  tr  ana  formation  matrix  to 
aatlafy  Eq  14. 

x'Xk)  * rV(k)  (14) 

The  control  law  synthesized  to  meet  Eq  (13)  is  as  follows. 

u(k)  = C'[h  + e]  T"x'(k)  (15) 

The  H and  L matrices  are  formed  from  the  Brimovsky  controllable 
canonical  form  of  AF  and  AG  respectively  (Ref  6:44-46).  The  Z 
matrix  Is  uaed  to  specify  the  desired  closed-loop  eigenvalues.  For 
tha  tlae-optlmal  control  law  the  closed-loop  eigenvalues  are  placed 
at  tha  origin  In  the  z-plane  and  the  E matrix  Is  the  zero  matrix. 

Infer  to  Reference  6,  page  45 • for  the  algorithm  uaed  for  calculating 
2 for  non-zero  eigenvalues. 

Matrix  Hanlpulatlon  Routines 

Table  II  shows  the  matrix  manipulation  routines  called  in  the 
original  FORTRAC  package  and  the  corresponding  routines  used  In  the 
version  edited  to  run  on  the  0)0  computers  at  Wrlght-Patterson  AFB. 
The  elganstructure  analysis  routines  were  taken  from  the  EISPACK 
Library  (lef  9)  available  on  the  CDC  computers.  The  other  routlnea 
mere  written  by  the  author  and  the  source  code  la  shown  In  Appendix  A. 
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Table  II 

Matrix  Manipulation  Routines 


Routines  Called 
In  The  Original 
rOSXRAC  Package 

FOICAF 

FOICBF 

FOICCF 

FOICDF 

FOICEF 

FOICJF 

F04AEF 


FOIATF 

FOIAKF 

F02AGF 

F02APF 


Routines  Used  In 
The  Wrlght-Patterson 
Adaptation 

PRESET 

IDEMT 

COPTAB 

ADDMAT 

SUBMAT 

TRANPOS 

INVERT 


Elgens tructure 
Routines 


BALANC 

KTJmRS 

ELTRAN 

HQR 

BQR2 

BALBAE 


Appendix  C gives  the  syntax  required  to  use  the  author's  matrix 
■anlpulatlon  routines  and  gives  a brief  description  of  the  algorithms 


used. 

The  above  routines  along  with  OISCLS  are  Included  In  the  INTERAC 
package.  The  subroutine  DISCLS  calculates  the  closed-loop  equations 
which  can  then  be  used  In  a discrete  simulation. 


L 


III.  Specialized  INTERAC  Software 

Before  adapting  FORTRAC  to  an  Interactive  program,  several 
Interactive  programs  (Ref  1, 2,3,4,  and  5)  available  on  the  Wrlght- 
Petterson  APB  computer  were  examined  to  see  how  efficient  they  were 
for  the  user.  The  first  Item  of  note  was  that  each  program  had 
different  Input  format,  output  content  and  format,  and  control 
structure.  This  greatly  Increases  the  workload  on  the  user  who  has 
to  use  several  programs  In  a design  effort.  With  the  exception  of 
OPTCGN,  there  are  no  provisions  for  correcting  an  erroneous  data 
point  entry  without  reentering  the  entire  data  set.  All  programs  are 
terminated  by  the  SCOPE  operating  system  for  a simple  Input  error  by 
the  user.  All  programs  allowed  the  user  to  select  the  options  he 
wished  to  execute,  but  there  are  no  provisions  to  stop  between  options 
and  allow  the  user  to  decide  If  the  sequence  should  be  continued.  In 
general  each  program  has  some  good  and  many  bad  points  In  regard  to 
user  efficiency. 

Three  methods  of  developing  PORTRAC  Into  an  efficient,  user 
oriented.  Interactive,  package  were  considered.  The  first  was  to 
redefine  the  operating  system  for  Interactive  computer  use.  Although 
this  would  provide  the  greatest  Improvement  for  all  classes  of 
Interectlve  programs,  the  time  and  resources  required  would  be 
beyond  the  capability  of  this  Investigation.  The  second  method  would 
be  to  utilise  a higher  order  computer  language  other  than  FORTRAN. 

Of  the  languages  avellsble  on  the  Wrlght-Petterson  CDC-6600/Cyber-74 
computer  system,  POBTRAN  Extended  seamed  to  be  the  beat  suited  for 
the  job.  This  left  direct  POBTRAN  prograemlng,  possibly  sacrificing 
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a Uttla  coaputar  rasource  afflclency  In  ordar  to  achlava  tha  dasirad 
usar  afflclancy. 


Tha  IHTESAC  softwara  dascrlbad  In  this  chaptar  was  davalopad  to 
Mat  tha  spaclflcatlons  glvan  In  Chaptar  1 using  tha  mair-ttninw  . 
capabllltlas  of  tha  FORTRAN  Extended  language  (Raf  10)  and  tha  SCOPE 
operating  systaa  (Ref  11)  available  on  tha  Wright-Pattarson  AFB 
eoaputars.  Tha  three  areas  in  a typical  program  operation  cycle 
addressed  are:  (1)  input,  (2)  output,  and  (3)  program  sequence 
control. 

Input  Includes  both  numeric  data  and  alphanumeric  commands.  The 
two  subroutines  davalopad  for  this  are  READNUM  (read  numeric)  and 
SBADCCM  (read  command) . The  data  output  is  controlled  in  content 
by  an  entry  in  a variable  output  table  and  in  format  by  the  subroutines 
PRIBTR  and  LISTER.  The  program  sequence  control  is  based  on  entries 
In  a valid  command  table  and  is  achlavad  in  this  program  with  tha 
subroutine  COMND. 


Input  Subroutines 

Tha  subroutine  READNUM  was  first  conceptualized  due  to  tha 
dafielancias  in  the  SCOPE  operating  systam  on  tha  Wright-Pattarson  AFB 
eoeputer.  Whan  a format  error  is  made  in  a numeric  input  operation,  an 
error  Mssage  is  printed  and  tha  program  is  terminated.  This  can  ba 


quite  frustrating  to  a user  who  is  Just  entering  the  last  number  of  a 


10  by  10  Mtrlx.  The  READNUM  subroutine  reads  all  Inputs  as  alphanumeric 
characters  snd  checks  for  forMt  errors.  If  any  errors  are  detected 
the  user  is  asked  to  reenter  those  iteM.  The  error  checking  is 
baaed  on  a lagal  entry  table  to  be  discussed  later. 


Tabic  III 

READNUMy  Legal  Entry  Table 


♦ft  ♦B 

♦C 

LEGAL  ENTRY  TABLE 

♦A 

♦B 

♦C 

LEGAL  ENTRY  TABLE 

♦♦  ♦♦ 

01 

66666666666666666666 

1 00 

02 

11111111111153555555 

5 

40 

34 

11111111111113111112 

ft  01 

03 

11111111111131555553 

6 

41 

35 

11111111111113111112 

B 02 

04 

11111111111151555555 

1 

42 

36 

11111111111113111112 

C 03 

05 

11111111111131555553 

S 

43 

37 

11111111111113111112 

D 04 

06 

11111111111151555555 

9 

44 

38 

11111111111113111112 

E 05 

07 

1111111111 1151555225 

+ 

45 

39 

11111111111113551552 

F 06 

08 

11111111111151555555 

- 

46 

40 

11111111111113551552 

6 07 

09 

11111111111151555555 

♦ 

47 

41 

11111111111113555552 

H 10 

10 

11111111111151555555 

50 

42 

11111111111142455554 

I 11 

11 

11111111111151555555 

< 

51 

43 

11111111111153555555 

J 12 

12 

11111111111151555555 

52 

44 

11111111111153555555 

K 13 

13 

11111111111151555555 

■1 

53 

45 

11111111111135777773 

L 14 

14 

11111111111131555553 

s= 

54 

46 

11111111111153555555 

M 15 

15 
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BBB^  I 
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START 


SITUATI(»I  1 


IS  ^ 
VARIABLE 
REAL  OR 
.DiTEGER^ 


LNTEGER 


CALL  SUBREL 


CHARGE  STORAGE  OF 
VARIABLE  VALUES  FROM 
INTEGER  TO  REAL 
REPRESENTATION 


IN(K) ,K-1,72 


72  CHARACTER  CARD 
IMACZ  BEAD  IN  R1 
FORMAT 


K-R-t-l 


ENTER  TABLE  WITH 
IN(K)-I>2  TO  FIND 
VALUE  FOR  GO 


Rl'i') 


SED)  7 


END  OP 

PRINT 

DECODE 

DECODE 

DECODE 

CARD 

ERROR 

INDICES 

ALPHANUMERIC 

NUMERIC 

DUGS 

AND  RESET 
POINTER 

COMMANDS 

ENTRIES 

flfim  2 


Flow  Chart,  Part  I 


flgun  3.  SEADNUM  Flow  Chart,  Part  II 


ntur*  5.  REAHnm  Flotf  Chart,  Fart  IV 
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PRINT  APPROPRIATE 
ERROR  MESSAGE 


^ IS^^ 
VARIABLE 
REAL  OR 
INTEGER 

'V  ? ^ 


INTEGER 


IS  ^ 
VARIABLE 
REAL  OR 
INTEGER^/ 

iniTEGER 


CALL  SUBINT 
CHANGE  STORAGE  OF 


NORMAL 


Flgur*  6.  READNtJM  Flow  Chart,  Part  V 
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Thtt  power  in  SEAmUM  is  particularly  evident  In  matrix  data  entry. 
Bach  element  of  the  matrix  Is  constantly  available  for  data  entry  or 
change.  A data  entry  error  noted  after  a carriage  return  can  still  be 
corrected  on  the  next  line.  Six  special  commands  give  the  user 
coaplete  control  over  the  data  Input  process. 

The  basic  flow  chart  for  REAMUM  Is  shown  In  Figure  2.  As  can 
be  seen,  this  subroutine  can  handle  variables  In  real  or  Integer 
representation.  All  data  within  the  subroutine  are  represented  as 
real  nusibers,  and  just  prior  to  termination  the  subroutine  converts 
thee  back  to  Integer  form  If  necessary. 

All  Input  characters  are  checked  against  the  legal  entry  table, 
Table  111.  This  table  Is  stored  as  an  array  of  65  computer  words. 

Bach  octal  bit  of  the  word  Is  used  In  a different  situation.  Thus 
each  character  can  cause  one  of  8 actions  In  up  to  20  different 
situations.  Situation  1 (the  table  Is  read  from  right  to  left)  Is  used 
at  tha  first  decision  point  shown  In  Figure  2.  The  resMlnder  of  the 
flow  chart  for  the  REAENOM  subroutine  Is  shown  In  Figures  3 thru  6. 

All  aatrlx  elements  are  entered  by  rows  In  free  format  until  all 
elements  have  been  filled,  unless  the  sequence  Is  altered  with  one  of 
the  special  commands.  The  first  special  character  of  Interest  is  the 
asterisk  (*),  shown  In  Figure  3.  The  entry  of  the  * causes  the  data 
entry  pointer  to  move  to  the  next  element  of  the  matrix  with  no 
ehanga  to  tha  valua  of  the  current  element.  Figure  4 shows  the  section 
for  decoding  the  six  special  alphanumarlc  conunds  and  the  special 
character  dollar  sign  ($).  Any  time  the  $ la  encountered  a return 
abort  Is  executed  from  the  subroutine. 

The  six  alphanumeric  commands  are  "abort**, "zero'*, "continue". 
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''nMd'',''ll8t'',  and  "/n,m/''.  The  "abort"  command  haa  the  same  result 
as  the  $.  "Zero"  causes  all  elements  of  the  variable  to  be  set  to 
aero  and  la  useful  for  entering  sparse  matrices.  "Continue"  results 
In  a normal  termination  of  the  subroutine  regardless  of  the  data 
entry  pointer  position.  "Read"  allows  the  Input  data  to  be  read  from 
the  file  TAFElf  TAFE2,  or  TAPES.  An  optional  parameter  following  the 
coHsand  "read  n"  Is  "rewind",  which  causes  TAFEn  to  be  rewound  prior 
to  reading  data  from  It.  The  command  "list"  catises  the  current  contents 
of  the  variable  elements  to  be  printed  at  the  terminal.  The  optional 
parameter  "/n,m/"  following  "list"  Is  used  to  print  only  one  element, 
one  row,  or  one  column  of  the  variable.  All  of  the  above  commands 
WKf  be  abbreviated  to  as  little  aa  the  first  letter.  The  command 
"/n,m/"  la  used  to  reset  the  data  entry  pointer  to  element  (n,m). 

This  allows  previous  errors  to  be  corrected  or  several  elements  that 
need  not  be  entered  to  be  skipped. 

The  technical  aspects  for  using  this  and  other  INTERAC  subroutines  ' 
for  different  programs  can  be  fotmd  In  the  Programmer's  Guide, 

Appandlr  C. 

The  second  Input  subroutine  developed  Is  READCOM.  This  siibroutlne 
{ reads  alphanumeric  commands  of  9 characters  or  leas  and  checks  them 

against  a list  of  valid  commands.  Commands  nay  have  any  characters, 
including  single  embedded  blanks,  except  the  slash  (/),  coomu  (,),  and 
dollar  algn  ($).  Any  abbreviation  that  la  unique  In  the  comaund  list 
' may  be  used.  The  flow  chart  for  READCOM  la  shown  In  Figure  7. 

Maltiple  commands  may  be  entered  If  they  are  separated  by  "/"  or  at 
least  tvo  blanks.  A command  string  la  also  recognised  and  Identified. 

A coanmnd  string  la  a coanmnd  followed  by  one  or  more  ocher  coamumda 
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acting  aa'nodlflng  pcranatars.  Tba  paraoateri 
coaBand  and  fron  each  other  by  conaws. 


iparatad  from  tha 


Tbara  ara  two  apaclal  coaauBida  that  ara  usad  intamally  in  tha 
anbrotttlna.  Tha  $ causes  a return  abort.  Anytlna  tha  $ Is  ancountarad 
tba  subroutine  returns  1th  sn  Indication  that  no  coinmands  have  been 
antarad.  Tha  second  spaclal  coansnd  is  "***••.  readccXI  normally  reads 
only  on^  80  character  input  line.  But  hen  ”***''  la  encountered  as 
a cflanwad  or  parasMtar,  another  line  is  read.  Thera  are  four  error 
conditions  recognised  by  READCOM.  After  the  error  message  is  printed 
tha  user  Is  asked  to  reenter  and  is  allowed  another  80  character  line. 

Output  Subroutines 

Five  of  the  legal  comnands  Interpreted  by  READCOM  are  used  to 
change  parameters  in  the  output  subroutines.  Three  of  these  commands 
ara  concerned  with  changing  entries  in  the  variable  output  table.  Data 
output  Is  sent  to  two  possible  locations  as  prescribed  by  the  user 
with  these  commands.  Data  can  be  displayed  at  the  terminal  or  can  be 

written  on  a local  file  (to  be  disposed  to  a line  printer  after  the  | 

program  is  terminated)  , or  both.  | 

Tha  subroutine  PRINTR  Is  called  by  the  coBq>utatlonal  subroutines  | 

whenever  the  value  of  a variable  is  changed.  In  this  subroutine  the  | 

variable  output  table  Is  checked  to  see  where  the  variable  values  \ 

should  be  printed.  If  display  at  the  teradnal  or  write  to  file  are  ^ 

Indicated  the  subroutine  LISTER  is  called.  Within  LISTER  the  data  is 
formated  to  appear  on  a 72  character  line  for  the  terminal  or  a 132 

I 

character  Una  for  the  Una  printer  output  as  appropriate.  The  number 
of  significant  digits  usad  in  each  caaa  can  be  specified  by  the  uaar. 
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Progran  Sequence  Control 


Tha  output  format  and  contant  as  well  as  all  program  operatlona 
ara  controllad  by  the  user's  alphanumeric  Inputs.  Table  IV  lists  all 
legal  conmands  available  to  the  user. 


Table  IV 

Legal  Commands  For  The  INTESAC  Package 


ntlllty 

Options 

I-O 

Options 

Design 

Options 

Variables 

OPTI(»IS 

ENTER 

RUN 

AMATRIX 

AFMATRIX 

CEIGEN 

INTEGRATE 

VASIABLES 

CHANGE 

FORTRAC 

BMATRIX 

AGMATRIX 

DEIGEN 

MODAL 

SAVE 

DFOBMAT 

SAMPLE 

RMATRIX 

AEMATRIX 

DESIGEN 

INVMODAL 

STQB 

OFOBMAT 

AUQIAT 

CMATRIX 

ARMATRIX 

CLEIGEN 

UREFORM 

EHD 

DISFLAI 

TRANSFORM 

FMATRIX 

ACMATRIX 

TSAMPLE 

APLFORM 

lESTART 

OUTPUT 

CONLAH 

OlATRIX 

BFMATRIX 

DIMENSION 

BRNFORM 

BEHIND 

DELETE 

CLOSELOOP 

DEMATRIX 

BQIATRIX 

STATES 

CINDICES 

PBINT 

SETDIMEN 

DRMATRIX 

TINVERSE 

CONTROLS 

COFMATRIX 

• 

KMATRIX 

COMMANDS 

CUIAIRIX 

DISTURBS 

Bafar  to  Appendix  B for  the  description  and  use  of  each  coonand.  The 
table  of  legal  commands  Is  stored  In  the  array  LIST,  dimensioned 
n by  7.  Col^a■n  1 of  LIST  contains  the  Hollerith  representation  of 
each  option  or  variable  name. 

The  other  columns  of  LIST  contain  all  the'  Information  needed  to 
decide  what  action  to  take  when  the  command  Is  Input  by  the  user. 
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Coltan  2 divides  the  commands  Into  four  types:  (1)  option  names, 

(2)  real  variable  names,  (3)  Integer  variable  names,  and  (4)  output 
variable  names.  The  output  variables  are  Intermediate  results  In  a 
computational  algorithm.  These  results  are  stored  In  temporary 
scratch  locations  in  memory,  and  the  user  may  not  enter  or  change 
data  within  them.  The  user  nay  Indicate  If  and  where  the  variable 
results  should  be  printed  when  they  are  calculated.  All  real  and 
Integer  variables  are  stored  In  contiguous  locations  In  a common 
block.  The  other  five  columns  of  LIST  Indicate  the  starting  location 
for  the  variable  within  the  common  block,  the  maximum  row  and  column 
dimensions  allowed,  and  the  currently  used  row  and  column  dimensions. 
Option  names  are  divided  Into  five  categories  by  an  entry  in  the 
last  column  of  LIST.  These  categories  are  shown  In  Table  V. 


Table  V 

Option  Catagorles 


latry  In  Column 

7 Of  LIST 

Description 

Principal  command  only,  with  n required 

parameters 

0 

Principal  command  only,  parameters  optional 

-1 

Principal  comnand  only,  no  parameters  allowed 

-2 

Principal  comnand  or  parameter,  no  parameters 
allowed 

-3 

Parameter  only 

Sobtoatlne  COMND  tiaes  the  entries  In  the  array  LIST  to  Initiate  the 
aetiom  requested  by  the  user's  comsand  Input.  Refer  to  Appendix  A for 
the  UstlAg  of  this  subroutine. 
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IV.  C-lAl  and  F-4E  Hathematlcal  Modeling 

An  analog  fly-by-wlre  control  system  for  the  C-141  has  been 
designed  and  built  by  Honeywell,  Inc.  under  contract  to  the  Air 
Force  Flight  Dynamics  Lab  (Ref  12) . The  results  of  this  project 
could  provide  data  against  which  to  check  a digital  control  system 
dealgned  using  INTERAC.  The  second  plant  chosen  to  test  the  INTERAC 
package  la  an  F-4E  (Ref  13) . This  chapter  describes  the  mathematical 
nodallng  process  for  the  two  plants.  After  several  computer  runs, 
the  C-141  was  dropped  from  consideration  due  to  modeling  problems 
Co  be  dlsaissed  later.  The  chapter  ends  with  a discussion  of  the 
raaults  obtained  for  the  F-4E. 

C-141  Model 

The  nodel  for  the  C-141  was  derived  using  data  from  Reference  12. 


The  longitudinal  equations  of  motion  are  given  as 

S = M^w  ♦MvvW  ♦MqG  + M/g  Xe +M;sp^sp  (16) 

w=  Zy  u Z^w  + Uo©  + Z/e  /e  Z/sp/sp  (17) 

u = Xy  u ♦ XvvW  - g©  + X/sp/sp  (18) 

h*  -w  ^ Uo©  (19) 

^Zpiiot  = w U)®  - 0 (20) 
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A A 

Th«  C aatrlx  was  than  chosan  to  giva  a C output.  C la  dascribad 


la  Baf  14  and  tha  aquation,  usad  la  (Rsf  12:55) 

' ^Zpitot  ^ (24) 

Kq  = 217  for  Uco=  350 


This  was  tha  first  hint  that  tha  data  provldad  by  Haf  12  nay  not  ba 
▼slid.  Using  tha  daflnltlon  of  C*  from  Raf  14,  tha  valua  of  Kq  sho\U.d 
haws  baan  10.8.  Tha  C natrlx  than  bacaaa. 


^ Is  442.2  or  452.8.  Ths  control  law  for  a discrete- time 
tlaa-optlmel  tracker  was  synthesised  for  both  C matrices  with 

4 

similar  results.  The  gains  were  excessive  (greater  than  10  ) and  a 
time  simulation  showed  huge  overshoots  and  Impossible  control 
requirements  (elevator  deflections  of  400*) . A modeling  problem 
was  suspected  when  It  was  noted  that  the  calculated  Brunovsky 
transformation  matrix  was  Incorrect.  This  was  due  to  computer 
numerical  problems  caused  by  the  fact  that  the  plant  as  modeled 
was  almost  uncontrollsble.  In  an  attempt  to  verify  the  modeling,  the 
coefficients  taken  from  Page  28  and  31  of  Ref  12  were  checked  against 
the  analog  flow  chart  na  Page  33  of  Ref  12.  There  were  numerous 

numerical  differences  and  a sign  difference  In  the  ( _ term  In  the 

•sp 

tt  equation  (Eq  18).  Thus  the  validity  of  any  comparison  of  a discrete 
fly-by-wlre  system  developed  using  INTERAC/FORTRAC  with  the  analog 
system  described  In  Ref  12  was  In  serious  doubt.  At  this  time  It 
was  decided  to  attestpt  a control  system  design  for  another  aircraft. 


f-4B  Model 

The  T-4E  was  modeled  using  data  from  a lab  project  from  the  EE  6.41 
course  at  A7IT  (Ref  13) . The  coefficients  of  Table  VII  were  tised  In  the 
following  equations  of  motion  t 
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Dlaanalonallzlng  and  other  ralatlonahlps 


• - w/g  C - -«g/Sq  u'-  u/H  « . 57.3*  Q - 57.30^ 
•'  The  elevator  actuator  and  aervo  wore  nodeled  as  the  first  order  lag 


Se  _ 20 

ejL  " s ^ 20 


(29) 
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Th«  £ matrix  was  chosen  to  allow  tracking  of  flight  path  angle,  @-e<, 
as  would  be  desirable  In  an  automatic  approach. 


c - [01-1000] 


Q2) 


.1 

I 


A control  law  was  synthesized  to  track  flight  path  angle  for  various 
sa^>llng  times . 

Bsaults  of  F-4E  Control  Law  Synthesis 

At  sll  sampling  times  a simulation  showed  the  F-4E  to  be  on 
track  within  4 sample  periods.  This  requires  higher  gains  and  higher 
control  power  as  the  sampling  Interval  decreases.  Table  VIII  shows 
some  of  the  representative  data  obtained  for  sampling  Intervals  of 
.05  to  5.0  seconds,  with  a unit  step  command  of  1 degree. 


Table  VIII 

Simulation  Data  for  the  F-4E  with  a Discrete-Time 
Flight  Path  Angle  Tracking  Control  System 


Haxlmum Final 

Saa^Ung  Elevator  Change  In  Negative  Magnitude  of  Change  In 
Interval  Deflection  Forward  Peak  In  Feedback  Forward 

Velocity  Commanded  Gain  Required  Velocity 
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At  saople  times  below  .5  seconds,  the  elevator  deflection  required  Is 
beyond  the  capability  of  the  aircraft.  A sample  time  of  2 seconds 
requires  the  least  control  power  and  produces  the  smoothest  response. 
But  8 seconds  to  track  a 1 degree  flight  path  change  Is  excessive  for 
aircraft  control.  Even  the  lowest  practical  sampling  time  of  .5 
seconds,  may  be  too  Ugh  for  this  fighter  type  aircraft. 

; I The  discrete-time  time-optimal  control  seems  to  work  quite  well , 

^ If  a relatively  long  sampling  Interval  can  be  tolerated.  If  a higher 

; j sai^llng  time  Is  required,  a sub-optimal  control  must  be  used.  One 

I possibility  Is  to  use  deadbeat  control,  which  In  this  case  would  allow 

7 sample  periods  to  track  the  command.  The  deadbeat  algorithm  Is  not 
presently  In  the  IMTESAC  package,  but  It  could  be  easily  added.  The 
second  possibility  Is  to  place  the  desired  eigenvalues  at  some  place 
other  than  the  origin  In  the  z-plane.  Since  there  are  an  Infinite 
number  of  possible  eigenvalues  to  use , the  INTESAC  package  Is  more  . 
desirable  than  the  batch  mode  program  for  quickly  trying  the  most 
promising  choices  of  eigenvalues. 
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V.  Conclusions  and  Reconmendatloos 

Conclusions 

Tha  INTERAC  package  meets  all  of  the  specifications  given  In 
Chapter  1.  There  Is  still  a little  room  for  Improvement  In  the  I-O 
versatility,  but  what  Is  now  presented  Is  a vast  Improvement  over 
existing  programs.  In  a current  Investigation,  Lt.  Stanley  J.  Larimer 
Is  Incorporating  many  of  these  user  oriented  Ideas  Into  a classical 
control  design  package. 

The  use  of  discrete-time  time-optimal  control  systems  for  total 
aircraft  control  Is  not  generally  possible  due  to  the  long  sampling 
Interval  required  to  keep  the  magnitude  of  the  control  within  practical 
limits.  There  are  applications,  with  low  frequency  disturbances  and 
comswnds,  where  the  time-optimal  control  could  be  used,  such  as  cnxlse 
control  and  possibly  automatic  approach  control.  The  total  aircraft 
control  system  will  probably  require  a sub-optimal  control  law. 

The  theory  behind  the  FORTRAC  and  INTERAC  packages  does  not 
provide  a cut  and  dry  control  system  design  process.  There  are  too 
many  possible  parameter  changes  In  the  multi-input,  multi-output 
system  for  a one  best  answer  to  be  achieved  In  one  Iteration.  The 
IBTESAC  package  provides  a basis  for  quickly  Investigating  the 
various  parsmster  changes. 

Escoisndatlons 

(1)  A graphical  simulation  capability  should  be  added  to  the 
URUUC  package.  The  magnitude  of  the  feedback  gains  does  give  a 
ganaral  idea  as  to  the  reaponce  of  the  closed  loop  system.  But,  as 
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•havn  for  th«  sampling  Intarval  of  .6  and. 7 saconds  in  Tabla  VIII, 
thmra  can  ba  a drastic  change  In  the  total  system  response  without 
large  changes  In  the  feedback  gains. 

(2)  The  LEADHEAJ  subroutine  for  achieving  the  sub-optimal  control 
law  should  be  added  to  the  package.  This  subroutine  Is  In  the  FOBTSAC 
package,  but  It  Is  not  called  by  the  master  program  furnished. 

(3)  The  investigation  of  the  discrete-time  fly-by-wire  control 
ayatem  for  the  C-141  should  be  continued  with  basic  flight  test  data 
from  Lockheed.  The  heavy  transport  type  aircraft  does  not  normally 
require  high  frequency  control  and  response,  and  thus  It  would  be 
Ideal  for  further  Investigating  the  time-optimal  control  laws. 
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INIERAC  Is  an  interactive  software  package  for  synthesizing  a 
discrete  state-variable  feedback  gain  matrix  to  control  a multi-input, 
multi-output  continuous  plant  described  by  Eq  (1) . 


x(t)  • A x(t)  + B_  u(t)  + R ^(t) 
j(t)  - C x(t) 


(1) 


where 


A is  an  N by  N state  matrix, 

7 is  an  N by  M control  matrix, 

R is  an  N by  HDD  disturbance  matrix,  and 
7 is  an  (NCC  or  NDD)  by  N output  matrix. 


y(t)  defines  the  NDD  outputs  that  are  to  reject 
disturbances  or  the  NCC  outputs  that  are  to  track 
input  commands. 


Three  types  of  problems  can  be  examined:  (1)  the  regulator 
problem I (2)  the  disturbance  rejector  problem,  and  (3)  the 
trecklng  problem.  The  computer  input  of  the  problem  type  and  the  plant 
deecriptlon  matrices  can  be  accoaq>llshed  directly  by  the  user  or  in 
response  to  requests  from  the  computer.  23  alphanumeric  options  are 
need  In  response  to  the  prompt  "COMMAND  »”  to  control  data  Input, 
decs  output,  and  program  aequence  control. 


Commands 

To  execute  the  program  for  the  first  time,  type  FORTRAC.  All 
further  input  data  will  be  requested  by  the  computer.  To  terminate 
the  program  type  STOP  In  response  to  the  prompt  "C(A1MAND»". 


i3 


s 1 

■f  J 


For  the  do-lt-yours«lf  user  the  following  comnands  are  available 
for  complete  Input, output,  and  program  sequence  control. 

(1)  raTER  [,  parameters]  This  command  allows  the  user  to 
•pecify  the  dimensions  of  a matrix  and  then  to  enter  the  values  for 
the  matrix  elements.  The  parameters  are  variable  names  and  their  use 
la  optional.  If  the  variable  names  are  not  Included  they  will  be 
requested  by  the  computer.  The  legal  variable  names  are  described 
under  options  16  through  21. 

(2)  STOP  This  connand  terminates  the  program  and  prints  a 
termination  message  Informing  the  user  of  the  contents  of  the  files 
AHSHER  and  DATA  that  are  left  by  the  program. 

(3)  EWD  This  command  has  two  functions.  First  In  response 

to  the  prompt  "COMMAND  the  program  Is  terminated,  but  the 

termination  message  Is  not  printed.  END  Is  also  used  to  terminate 
sub-options  such  as  when  the  command  ENTER  Is  given  without  a parameter 
list. 

(4)  PRINT  [,  parameters]  This  command  causes  the  current 
valnes  of  the  variables  Identified  In  the  parameter  list  to  be  printed 
at  the  terminal  and  to  be  written  to  the  file  ANSWER.  The  parameters 
are  optional , and  If  not  usad  the  variable  names  will  be  requested 

by  the  cos^uter. 

(5)  DPORMAT  This  command  allows  the  user  to  specify  the  number 
of  significant  digits  to  be  used  fox  numbers  printed  at  the  terminal. 

(6)  [,  paramaters]  This  conund  allows  the  user  to 


change  the  values  stored  In  a variable.  The  dimensions  of  matrices 
cannot  be  changed  with  this  command.  The  parameters  are  variable  names 
and  their  use  la  optional.  If  the  parameters  are  not  Included  the 
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n«—  will  be  requested  by  the  computer. 


(7)  DISPLAY  [,  parameters]  This  is  a dual  function  command. 
With  the  parameters  (variable  names) , the  variables  to  be  printed  at 
the  terminal  are  specified.  Without  the  parameters,  the  variables 
already  designated  to  be  printed  at  the  terminal  are  Identified  and 
'the  names  are  printed.  The  values  of  the  variables  on  the  display 
list  are  printed  at  the  terminal  whenever  the  variable  Is  first 
encountered  or  when  Its  value  Is  changed  In  the  control  design 
algorithms . 

(8)  OUTPUT  [,  parameters]  This  command  Is  used  the  same 

as  DISPLAY,  except  OUTPUT  refers  to  the  variables  to  be  written  to 
the  file  ANSWER. 

(9)  DELETE  [,  parameters]  This  comnand  deletes  the  named 

variables  (parameters)  from  both  the  display  and  output  lists.  If  no 
parameters  are  given,  then  all  variables  are  removed  from  both  lists. 

(10)  OPTIONS  This  command  causes  a list  of  the  valid  options 
to  be  printed  at  the  terminal. 

(11)  VARIABLES  This  command  causes  a list  of  the  legal 
variable  names  to  be  printed  at  the  terminal. 

(12)  SAVE  This  command  causes  the  current  contents  of  all 
variables,  all  matrix  dimensions,  and  the  display/output  list  to  be 
saved  on  the  file  DATA.  The  date  and  time  of  the  save  Is  printed,  and 
they  can  be  used  to  Identify  the  data  set  when  It  Is  reloaded  with 

I the  eoBoand  RESTART. 

(13)  RESTART  [,  REWIND]  This  command  causes  the  values 
saved  on  the  file  DATA  to  be  reloaded  Into  memory.  The  optional 
parasMter  REWIND  Is  used  to  rewind  the  file  before  reading  from  It. 
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Thtt  fll«  DATA  la  a sequential  file,  and  thus  RESTART  may  have  to  be 
used  several  times  to  get  the  desired  data  set  If  SAVE  was  used  more 
than  once.  The  date  and  time  that  the  data  set  was  saved  Is  printed 
to  Identify  the  set  that  Is  reloaded. 

(14)  REWIND  This  command  can  be  used  by  Itself  to  rewind  the 
file  DATA  or  as  a parameter  for  RESTART  as  shown  above. 

(15)  RUN  , parameter  This  command  causes  the  control  design 
algorithms  to  be  executed  starting  with  the  subroutine  named  In  the 
parameter  and  terminating  after  calculating  the  closed-loop  matrix. 

The  valid  parameters  are  FORTRAC,  SAMPLE,  AUGMAT,  TRANSFORM,  CONLAW , 
and  CLOSELOOP.  All  necessary  data  Is  assumed  to  have  been  entered 
and  all  matrices  are  assumed  properly  dimensioned. 

(16)  FORTRAC  This  command  by  Itself  causes  all  necessary 
data  to  be  requested  and  then  all  control  design  algorithms  to  be 
executed.  When  this  command  Is  used  as  a parameter  for  RUN,  all 
control  design  algorithms  are  executed  with  the  data  currently 

In  memory. 

(17)  SAMPLE  This  command  causes  the  sampled-data  system 
description  to  be  computed  from  the  continuous  system  description. 

The  variables  mxist  be  entered  prior  to  Issuing  this  command:  AMATRIX  - 
the  state  matrix,  BMATRIX  - the  control  matrix,  DIMENSION  - the  array 
containing  the  plant  dimensions  (N,  M,  NCC,  NDD),  RMATRIX  - the 
disturbance  matrix  only  If  NDD  la  non-zero  , and  TSAMPLE  - the 
desired  sampling  tlma  Interval.  Values  for  the  following  variables 
axe  calculated:  CEIGEN  - the  continuous  eigenvalues;  MODAL  - the 
modal  matrix;  INVMODAL  - the  Inverse  modal  matrix;  and  tha  dlscreta 
system  matrices  FMATRIX  - the  state  matrix , GMATRIX  - the  control 
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iutxix,  DBMATRIX  - the  disturbance  matrix  If  NBD  Is  non-zero , and 
DEIOt  - the  z - plane  eigenvalues . 
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(18)  AUGMAT  This  command  causes  the  sampled-data  system 

description  to  be  augmented  with  discrete  Integrators  for  disturbance 
rejection  or  tracking  problems.  These  variables  must  be  available 
In  memory  prior  to  Issuing' this  command:  FMATRIX,  GMATRIX,  DRMATRIX  - 
if  HDD  Is  non-zero,  CMATRIX  - the  output  matrix  If  NSD  or  NCC  Is 
non-zero,  DIMENSION,  and  INTEGRATE  - an  array  containing  the  number 
of  Integrators  desired  for  each  NDD  or  NCC  output.  One  Integrator 
is  used  to  track  or  reject  a step  Input,  two  Integrators  for  a ramp,  | 

etc.  Values  for  the  following  variables  are  calculated:  AEMATRIX-  * 

the  augmented  state  matrix,  AGMATRIX  - the  augmented  control  matrix,  ; 

1 

AIMAIRIX  - the  augmented  Input  command  matrix , ARMATRIX  - the 

augmented  disturbance  matrix,  and  ACMATRIX  - the  augmented  output  I 

matrix. 

(19)  TRANSFORM  This  command  causes  the  ^ and  ^ matrices  j 

from  the  augmented  system  description  to  be  transformed  Into  the  ' 

Brunovsky  controllable  canonical  form.  AFMATRIX,  AGMATRIX,  and  I 

DIMENSION  must  be  available  In  memory  prior  to  Issuing  this  command. 

Values  for  the  following  variables  are  calculated  : UREFORM  - the 
upper  row  echelon  form  of  the  matrix  [AG  , APLFORM  - the  Aplevlch 
form  .0  f the  matrix  [AG  AF] , BRNFORM  - the  Brunovsky  form  of  the 
matrix  [AG  AFl , BFMATRIX  and  BQIATRIX  - the  Bnmovsky  controllable 
canonical  form  of  AFMATRIX  and  AQIATRIX  respectively,  TINVERSE  - 
the  Inverse  transformation  matrix,  and  CINDICES  - the  controllability 
indices. 

(20)  CONLAW  This  command  causes  the  state-variable  feedback 


gain  matrix  K to  be  calculated  for  the  desired  closed-loop  eigenvalues. 
BFKAIRlXp  BGMATRIX,  TINVERSE,  DIMENSION,  and  DESEIGEN  - the  desired 
eigenvalues,  must  be  available  In  memory  prior  to  Issuing  this  command. 
Values  for  the  following  variables  are  calculated  : COFMATiHX  - the 
coefficient  matrix  which  will  be  the  zero  matrix  unless  non-zero 
desired  eigenvalues  are  entered,  and  KMA.TRIX  - the  feedback  gain 
matrix. 


(21)  CLOSELOOP  This  command  Is  used  to  calculate  the  closed- 
loop  state  matrix  and  the  closed-loop  eigenvalues.  The  variables 
required  In  memory  prior  to  Issuing  this  command  are:  AFMATRIX, 
AGMATRIX,  DIMENSION,  and  KMATRIX.  CLMATRIX  - the  closed-loop  matrix 
and  CLEIGEN  - the  closed-loop  eigenvalues  are  calculated. 

(22)  OFORMAT  This  command  allows  the  user  to  specify  the 
number  of  significant  digits  to  be  used  for  numbers  written  to  the 
file  ANSWER. 

(23)  SETDIMEN  This  command  properly  dimensions  all  matrices 
after  th.e  plant  dimensions.  In  the  array  DIMENSION,  have  been  entered. 

(24)  "variable  name”  When  the  variable  name  Is  typed  for 

a command,  the  user  Is  allowed  to  set  the  dimensions  of  the  variable 
and  to  enter  or  change  the  values  of  the  elements  of  the  variable. 

(25)  DIMENSION  This  command  opens  the  array  DIMENSION  and 
allows  the  user  to  enter  N - the  number  of  states,  M - the  number  of 
controls,  and  NCC  - the  number  of  commands  or  NDD  - the  number  of 
disturbances.  The  elements  of  DIMENSION  are  individually  addressable 
hj  the  variable- names  STATES,  CONTROLS,  COMMANDS,  and  DISTURBS. 

The  above  commands  are  used  In  response  to  the  prompt  "COMMAND  » " . 
If  tbs  command  has  paramaters,  the  parameters  must  be  separated  by 
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coaoMS.  Multiple  commands  may  be  entered  on  one  line  If  they  are 
separated  by  slashes  (/).  Execution  of  entered  commands  begins  after 
the  carriage  return,  but  to  delay  execution  and  enter  more  commands 
or  parameters, the  last  entry  on  the  line  must  be  "***''. 

When  a variable  is  opened  for  entering  or  changing  data,  there 
are  six  special  commands  available  to  facilitate  the  process.  The 
first  is  ABORT  which  will  terminate  current  commands  and  return  the 
prompt  "COMMAND »"  to  the  user.  ZERO  is  used  to  set  all  elements 
of  the  variable  to  zero  and  la  xiseful  for  entering  sparse  matrices. 

READ  n allows  the  user  to  Input  the  data  from  the  file  TAFEl,  TAPE2, 
or  TAPES.  The  optional  parameter  REWIND  following  READ  n,  is  used  to 
rewind  the  file  before  reading  from  It.  LIST  Is  used  to  list  the 
current  contents  of  the  variable.  The  optional  parameter  "/n,m/''  after 
LIST  is  used  to  print  only  one  element,  one  row,  or  one  column  of 
the  variable.  CONTINUE  closes  the  variable  and  starts  execution 
of  any  remaining  commands  or  Issues  the  prompt  "COMMAND»".  The 
last  special  command  Is  "/n,m/",  which  is  used  to  reset  the  data 
entry  pointer  to  element  (n,m).  All  matrix  elements  are  entered 
by  rows  In  free  format,  and  this  last  command  can  be  used  to  go 
back  and  correct  an  erroneous  entry  or  to  skip  over  elements  that 
need  not  be  entered. 

A dollar  sign  ($)  enterad  at  any  tlma  will  terminate  current 
activity,  cause  pending  comnands  to  be  Ignored,  and  will  ratum  tha 
, prompt  "COItfAND»"  to  the  user. 

The  remainder  of  this  guide  will  be  devoted  to  showing  examples 
of  the  use  of  the  sbove  commands. 


COnMAND  >> ENTER 


ENTER  VARIABLE  NAME  DR  END  »AMATRIX 


ENTER  DESIRED  DIMENSIONS  OF  AMATRIXf  MAX  ALLOWED  IS  < 2 Of  20> 


AHATRIX  DIMENSIONED  < 4*  4>  IS  OPEN 


LIST 


ANATRIX 


ANATRIX 


CONTINUE 


ENTER  VARIABLE  NAME  OR  END  »END 


COMMfiND  » ENTER »flMflTR IX 


ENTER  DESIRED  DIMENSIONS  OF  flMRTRIXf  MfiX  fiLLOUED  IS  < 20» 


flMRTRIX  DIMENSIONED  < 2»  2>  IS  OPEN 


COMMftND  > > PR I NT > RMfiTR I X/CHRNGE > RMRTR I X/PR I NT » RMRTR IX 


RNRTRIX 

l.OOE-^00 

2.00E^00 


2.00E-^00 

S.OOE^OO 


RNRTRIX  DIMENSIONED  < 2*  2>  IS  OPEN 

6>7*  8»  9 


RNRTRIX 

6.00E4>00 

S.OOE-^OO 


7.00E'»‘00 

9.00E-^00 


91 


I 


a ■ 


CDMHflHD  >>SflVE 

DftTR  SftVED  RT  12/13/77  17.15.37. 

CORMRNr  > > ENTER  > RMRTR I X^-RUN  » FDRTRRC 

ENTER  DESIRED  DIMENSIONS  OF  RMflTRIX»  MAX  ALLOWED  IS  < 20.  20). 


RMRTRIX  DIMENSIONED  < 3.  3)  IS  OPEN 

I.2.3.  0. 1»^3»2>'»2»L 


13  1.0000000 

23  0. 

33  4.0000000 


RMRTRIX 

2.0000000 

1.0000000 

2.0000000 


3.0000000 


2.0000000 


COMMAND  >> RESTART. REWIND 


DATA  SAVED  AT  12/13/77  17.15.37.  RELOADED. 


COMMAND  »PRINT.AM 


RMRTRIX 

6.00E4>00 

S.OOEfOO 


7.00E^00 

9.00E4>00 


COnMAND  » 01 SPLAY 


VARIABLES  TO  BE  DISPLAYED  AT  TERMINAL 
DEI6EN 
FMATRIX 
GMATRIX 
CLEI6EN 
KMATRIX 
TSAMPLE 
CINDICES 
BRNFDRM 


COMMAMD  >>DISPLAY» AMATRIX»BMATRIX 


COMMAND  »D1  SPLAY 

VARIABLES  TO  BE  DISPLAYED  AT  TERMINAL 
AMATRIX 
BMATRIX 
DEI6EN 
FMATRIX 
6MATRIX 
CLEI6EN 
KMATRIX 
TSAMPLE 
CINDICES 
BRNFORM 


COMMAND  » DELETE » TSAMPLE » CINDICES 


COMMAND  »DISPLAY 

VARIABLES  TO  BE  DISPLAYED  AT  TERMINAL 
AMATRIX 
BMATRIX 
DEI6EN 
FMATRIX 
6MATRIX 
CLEI6EN 
KMATRIX 
BRNFORM 
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COMMAND  » OUTPUT 


VARIABLES  TO  BE  PRINTED  TO  ANSWER 
AMATRIX 
6MATRIX 
DIMENSION 
RMATRIX 
CMATRIX 
FMATRIX 
GMATRIX 
AFMATRIX 
A6MATRIX 
AEMATRIX 
ARMATRIX 
ACMATRIX 
CLE16EN 
KHATRIX 
OESEIGEN 
MODAL 


COMMAND  >> DELETE 

ALL  VARIABLES  HAVE  BEEN  REMOVED  FROM  THE  DISPLAY  AND  OUTPUT  LISTS. 


COMMAND  >>OUTPUT^DISPLAY 


VARIABLES  TO  BE  PRINTED  TO  ANSWER 
VARIABLES  TO  BE  DISPLAYED  AT  TERMINAL 


o 


( ■) 


COMMAND  >> OPT IONS 
VALID  COMMANDS 


ENTER 

STOP 

END 

C» PARAMETER 

LIST] 

PRINT 

DFORMAT 

C» PARAMETER 

LIST] 

CHANGE 

C» PARAMETER 

LIST] 

DISPLAY 

C> PARAMETER 

LIST] 

OUTPUT 

C» PARAMETER 

LIST] 

DELETE 

OPTIONS 

VARIABLES 

SAVE 

C» PARAMETER 

LIST] 

RESTART 

REWIND 

C» PARAMETER 

LIST] 

RUN 

FORTRAC 

SAMPLE 

AU6MAT 

TRANSFORM 

CONLAU 

CLOSELOOP 

OFORMAT 

SETDIMEN 

» PARAMETER  <S> 

COMMAND  » 
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COMMAND  >> VARIABLES 

VALID  VARIABLES 
AMATRIX 
BMATRIX 
DEI6EN 
DIMENSION 

♦ CDFMATRIX 
INTEGRATE 
RMATRIX 
CMATRIX 
FMATRIX 
6MATRIX 
DEMATRIX 
DRMATRIX 
AFMATRIX 
A6MATRIX 
AEMATRIX 
ARMATRIX 
ACMATRIX 
CEIGEN 
CLEI6EN 
BFMATRIX 
BGMATRIX 
T INVERSE 
KMATRIX 
DESEIGEN 
CLMATRIX 
TSAMPLE 
STATES 
CONTROLS 
COMMANDS 
DISTURBS 

♦ MODAL 

♦ INVMODAL 

♦ CINDICES 

♦ UREFORM 

♦ APLFORM 

♦ BRNFORM 


MAX  DIMENSION  ALLOWED 


< 

20> 

20> 

< 

20* 

10> 

< 

20* 

2> 

< 

1* 

1> 

< 

0* 

0> 

< 

1* 

10> 

< 

20* 

5> 

< 

10* 

20> 

< 

20* 

20> 

< 

20* 

10> 

< 

20* 

5> 

< 

20* 

5> 

< 

20* 

20> 

< 

20* 

10> 

< 

20* 

5> 

< 

20* 

5> 

< 

10* 

20> 

< 

20* 

2> 

< 

20* 

2> 

< 

20* 

20> 

< 

20* 

10> 

< 

20* 

20> 

< 

10* 

20> 

< 

20* 

2> 

< 

20* 

20> 

< 

1* 

1> 

< 

1* 

1> 

< 

1* 

1> 

< 

1* 

1> 

< 

1* 

1> 

< 

0* 

0> 

< 

0* 

0> 

< 

0* 

0> 

< 

0* 

0> 

< 

0* 

0> 

< 

0* 

0> 

♦ ~ TEMPORARY  VARIABLE  PRINTED  ONLY  DURING  PROGRAM  EXECUTION. 


COMMAND  >> 


RMftTRIX 


13 

1.00000E-i>00 

2.00000E'«'00 

3.00000E+00 

4.  00000E-»>00 

21 

0. 

l.OOOOOE-^OO 

0. 

0. 

31 

4.00000E'»-00 

3.00000E-»-00 

2.  OOOOOE'i'OO 

l.OOOOOE+00 

43 

0. 

0. 

0. 

S.OOOOOE-^OO 

i 


I 

i ! 


I-  i 


BHaTRIX 


13 

S.OOOOOE-t-OO 

0. 

21 

0. 

0. 

33 

l.OOOOOE-i-00 

0. 

41 

0. 

1.00000E-i>00 

COMMflMD  >>DFaRriflT 

ENTER  NUMBER  QF  SIGNIFICANT  DIGITS  0 < N < 15. 


COMMAND  » PRINT 

ENTER  VARIABLE  NAME  DR  END  »AMATRIX 


ANRTRIX 


13 

1.00E-*>00 

2.00E'^00 

3.  00E-»>00 

4.00E'*>00 

23 

0. 

l.OOEfOO 

0. 

0. 

31 

4.00E‘*>00 

3.00E-»>00 

2.002-^00 

l.OOE+00 

43 

0. 

0. 

0. 

5.00E-i>00 

ENTER  VARIABLE  NAME  OR  END  »BMATRIX 


MATRIX 

t3  3.00E>»00  0. 

no.  0.  . 

31  . l.OOE^OO  0. 

43  . 0.  1.00E4>00 
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Thtt  following  pages  give  syntax  for  calling  the  specialized 
mTERAC  subroutines.  All  elements  in  the  call  statement  are  defined 
and  any  required  common  blocks  are  identified. 
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ROW  DIMENSION  OF  A AND  R 
COLUMN  DIMENSION  OF  A AND 


NOn  — ROM  OIHFNSION  OF  0 FROM  OIHFNSION  STATEMENT  IN 
CALLING  PROGRAM 


CALL  INVFRT(A,NA0,B,N9n,N,IFAlL)  «*  A « B (INVERSE) 

A -r  OUTPUT  MATRIX 

NAO  — ROW  OIHFN!?ION  OF  A FROM  DIMENSION  STATEMENT  IN 
CALLING  PROGRAM 


Th«  use  of  the  batch  mode  program  FORTRAC  Is  described  In  the 
preliminary  report  "FORTRAC:  A Software  Package  for  the  Design  of 
Multivariable  Digital  Control  Systems"  which  Is  available  from  the 
FGL  branch  of  the  Flight  Dynamics  Lab>  The  one  point  that  needs  to 
be  expanded  on  Is  the  card  sequence  required  for  the  different  types  of 
problems  that  can  be  run  with  the  program.  The  following  pages 
give  the  card  sequence  required  for  the  three  types  of  problems; 
the  regulator,  tracker,  and  disturbance  rejector. 

The  following  limits  are  In  effect  with  the  presently  dimensioned 
master  program. 

(1)  The  maximum  number  of  states  Including  augmenting 
Integrators  allowed  Is  7. 

(2)  The  maximum  number  of  controls  Is  4. 

(3)  The  maxlsaim  nuinber  of  augmenting  Integrators  allowed 
Is  3. 

(4)  The  maximum  nusd>er  of  samples  that  can  be  displayed 
In  either  the  discrete  or  continuous  simulation  Is  33. 
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Car4  S«qttMic«  for  the  Regulator 

1 ■•M JICC .HDD  iUYT  fiP? 

I a ntiab«r  of  states 
M • nuBiber  of  controls 
■CC  - 0 

no  - 0 

nr  " nuabex  of  outputs 
■PI*  - 0 


The  elements  of  the  state  matrix,  one  row  per  card 


The  elements  of  the  control  matrix,  one  row  per  card 


The  elesients  of  the  output  matrix,  one  row  per  card 


T The  sampling  Interval  to  be  uaed 


The  desired  closed-loop  eigenvalues , N cards  with  a real 
and  an  Imaginary  part  on  each  card 


0 No  observer  will  be  designed 


K An  Integer  Indicating  the  power  to  which  the  closed-loop 
matrix  shoxtld  be  raised 


NDS  Humber  of  time  Intervals  required  In  the  discrete-time 
simulation 


H Initial  conditions  for  the  states,  all  on  one  card 


1P8,H1S,NSII» 

Dg  • the  accuracy  to  be  used  In  the  Kutta-Merson  subroutine 
(suggest  approximately  .0000001) 

■D8  • nu^er  of  sample-time  Intervals  required  in  the 
continuous -time  simulation 

nZDS  ■ number  of  divisions  required  In  each  sample -time 
Interval 


■ conditions  for  the  states,  all  on  one  card 


i 

4 

I 
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Card  S«qu«nc«  for  the  TraclcT 


1 n^^ccjsnoijsnjst? 

I • nuid>«r  of  stmtes 
M ■ nuaber  of  controls 
HOC  " number  of  coomsnds 
HDD  - 0 

HIT  m number  of  outputs 

IPP  ■ nuDd)er  of  outputs  integrated  (equal  to  NCC) 


The  elements  of  the  state  matrix,  one  row  per  card 
The  elements  of  the  control  matrix,  one  row  per  card 


The  elements  of  the  connand  matrix,  one  row  per  card 
(this  must  be  entered  even  If  the  matrix  Is  null) 


The  elvnents  of  the  output  matrix,  one  row  per  card 


The  elements  of  the  matrix  defining  the  outputs  that  will 
track  the  Input  command 


T The  sampling  Interval  to  be  used 


The  nu^er  of  Integrators  desired  for  each  integrated 
output,  all  entered  on  one  card 


1 1 

i 


The  desired  closed*loop  eigenvalues , (M  -f  the  sum  of  the 
numbers  entered  on  the  preceeding  card)  cards  with  a real 
end  an  Imaginary  part  on  each  card 


0 Bo  (d>server  will  be  deslnged 


K da  Integer  Indicating  the  power  to  idilch  the  closed-loop 
matrix  should  be  raised 


IDS  lumber  of  time  Intervals  required  In  the  discrete-time 
slmuletlon 


) 


1 

8 

I 

] 

8 

! 


1 
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M initial  condltlona  for  the  states,  all  on  one  card 


Initial  conditions  for  the  Integrators  , all  on  one  card 


EPS,HDS,NSIDS 


Sane  as  for  the  regulator 


N Ixiltlal  conditions  for  the  states,  all  on  one  card 


Initial  conditions  for  the  Integrators , all  on  one  card 


Card  Sequence  for  the  Disturbance  Rejector 


1 N,M,MCC,NDD,NYY,NPP 

H ■ number  of  states 
M • number  of  controls 
HCC  - 0 

HDD  > number  of  disturbances 
NTY  ■ number  of  outputs 

HPP  m number  of  outputs  Integrated  (number  of  outputs  to 
reject  disturbances) 


The  elements  of  the  state  matrix,  one  row  per  card 


The  elements  of  the  control  matrix,  one  row  per  card 


The  elements  of  the  disturbance  matrix,  one  row  per  card 


The  elements  of  the  output  matrix,  one  row  per  card 


The  elements  of  the  matrix  defining  the  outputs  that  will 
reject  disturbances 


T The  sampling  Interval  to  be  used 


All  remaining  cards  are  the  same  as  for  the  tracker 
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